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Abstract. The flux from distant type la supernovae (SN) is likely to be amplified or 
de-amplified by gravitational lcnsing due to matter distributions along the linc-of-sight. 
A gravitationally lensed SN would appear brighter or fainter than the average SN at 
a particular redshift. We estimate the magnification of 26 SNe in the GOODS fields 
and search for a correlation with the residual magnitudes of the SNe. The residual 
magnitude, i.e. the difference between observed and average magnitude predicted by 
the "concordance model" of the Universe, indicates the deviation in flux from the 
average SN. The linear correlation coefficient for this sample is r = 0.29. For a similar, 
but uncorrclated sample, the probability of obtaining a correlation coefficient equal 
to or higher than this value is ~ 10%, i.e. a tentative detection of lensing at ~ 90% 
confidence level. Although the evidence for a correlation is weak, our result is in 
accordance with what could be expected given the small size of the sample. 
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1. Introduction 



Matter inhomogeneities in the Universe, such as galaxies, groups, and clusters of 
galaxies, give rise to gravitational fields which affect the path of light rays emitted 
from distant sources. Gravitational bending of light, so called gravitational lensing, can 
amplify or de-amplify the flux from distant sources. A gravitationally lensed standard 
candle, such as a type la supernovae (SN), at a particular redshift would thus appear 
brighter or fainter to an observer than the average standard candle at this redshift. 
Using observations of the matter in the foreground, we can estimate the magnification 
factor /i, which describes the amplification (fi > 1) or de- amplification (// < 1) relative 
to a homogeneous universe, of a line-of-sight. The amplified flux of a point source with 
flux / is given by fif. If our estimates of the amplification are correct, we expect to find 
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a correlation between standard candle brightness and magnification. In this letter we 
use SNe to search for this correlation. Gravitational lensing adds extra scatter, which is 
expected to grow with redshift [lj, to the intrinsic dispersion of SN brightness. However, 
most of this dispersion is due to a few highly magnified SNe and a large portion of the 
SNe are slightly de-magnified. The effects of intrinsic scatter and measurement errors 
are thus larger than the effect of gravitational lensing for most SNe and the correlation 
might consequently be hard to reveal. 

The correlation between foreground matter and SN brightness has been searched for 
earlier. Williams and Song [2] reported a correlation between luminosity of SNe and the 
density of foreground galaxies within 5 — 15 arcminutes radius from the position of the 
SN. Menard and Dalai [3] searched for a correlation on scales 1 — 10 arcminutes, but could 
not find any correlation. In this letter we search for a correlation between SN brightness 
and the magnification due to foreground galaxies within ~ 1 arcminute. Since we use 
spectroscopic or photometric redshifts of the galaxies, we take the spatial distribution 
of matter into account. The possible confusion of foreground and background galaxies 
have been incorporated in our error estimation. 

2. Data 

In order to reconstruct the amplification along a line-of-sight we need galaxy positions, 
redshifts, and luminosities. Hubble Deep Field North (HDFN) and Chandra Deep Field 
South (CDFS) observed within the Great Observatories Origins Deep Survey (GOODS) 
are ideal for our purposes, since they have been extensively imaged in many wavelength 
bands [HE] and have been extensively searched for SNe [6J, [Tl [8] . We use the information 
from the GOODS to estimate the magnification of 26 SNe observed in these fields. Since 
most of the galaxies lack spectroscopic redshifts (especially high redshift galaxies) we 
have to rely on photometric redshifts which depend on observations in many different 
wavelength bands. The galaxy data we have used have been described in detail in 
Section 3 in [9]. 

3. Estimating the magnification factor 

To compute the magnification factor of the SNe we have used the prescription described 
in [9]. Since each galaxy along the line-of-sight can influence the path of a light ray, we 
use a numerical code (QLET)[10] utilizing the multiple lens plane formalism and include 
all foreground galaxies within 1 arcminute from the SN in the calculations. Galaxies are 
modelled as truncated Singular Isothermal Spheres (SIS) or truncated Navarro, Frenk, 
and White profiles (NFW)[H]. The mass of the galaxy halos were estimated using 
the luminosities of the galaxies via empirical Tully-Fisher and Faber- Jackson relations. 
The discrepancy between the total mass of the galaxies in the fields and the total 
amount of matter (Qm ~ 0.3) in the Universe, is distributed homogeneously. Errors 
in the magnification factor of each SN were estimated using Monte-Carlo simulations, 
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Table 1. Residual magnitudes and magnification of supernovae. 
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Quoted errors are root mean square deviations from the mean. 



taking into account photometric redshift uncertainties and scatter in the Tully-Fisher 
and Faber- Jackson relations, following the method outlined in Section 4 in [9] . However, 
the largest uncertainty is due to the uncertainty in the halo model, which is reflected by 
the differences in fi calculated using SIS or NFW halo profiles. We have estimated the 
magnification factor for 26 SNe in the GOODS fields for which the distance modulus 
have been accurately measured^ [8]. For 14 of these SNe magnification factors were 
presented already in [9]. Magnification factors in logarithmic units, A /t = — 2.51og 10 /i, 
estimated for each SN assuming galaxy halos to be described by NFW- or SIS-profiles 
are listed in Table CD 

t We have only considered "gold" SNe. 
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4. The correlation 

Since flux is conserved, the average flux from a large number of lensed standard candles, 
all at the same distance from the observer, is expected to be the same as the flux from 
a standard candle in an homogeneous universe. The average brightness or distance 
modulus of a standard candle can thus be predicted from a cosmological model of 
a homogeneous universe. We take the "concordance model", a flat universe with a 
cosmological constant and matter density flu = 0.29, to be our model of the Universe. 
In Table [1] residual magnitudes, A, of the SNe calculated by subtracting the predicted 
distance moduli from the observed value are presented. SNe brighter than average thus 
have negative residuals, while fainter than average SNe have positive residuals. 

Figure [1] shows a scatter plot of residuals vs magnifications, the latter computed 
using NFW-profiles and expressed in logarithmic units A M = — 2.51og 10 /z. For the 
sample considered here the linear correlation coefficient is r = 0.29 (the correlation 
coefficient is slightly larger, r = 0.34, if the magnification is computed using SIS- 
profiles). If we instead compute a rank-order correlation coefficient (Spearman), we 
get a slightly lower value r s = 0.27 (r s = 0.31 for the SIS case). For a linear correlation 
between magnification and residuals, the data should give a reasonable fit to a straight 
line, A = a + bA^. Errors are large in both residuals and magnifications, therefore we 
take into account errors in both coordinated The best fit to the data is indicated by 
the solid line in Figure d For this line (a = -0.098 ± 0.048 and b = 1.4 ± 1.4) the 
chi-square is x 2 = 24.09 for 24 degrees of freedom. If there is no correlation we expect 
a line with a = and 6 = 0, indicated by the dashed line in the figure, to be a good 
fit to the data. The chi-square for this line is x 2 = 29.24, indicating a weak preference 
(~ 90% confidence level) for a positive correlation. 

5. Discussion 

The evidence for a correlation between residuals and magnifications is rather weak. 
Considering the low statistics, what could we expect to measure? Even if there 
exists a complete positive correlation between magnifications and residuals, the intrinsic 
dispersion in SN luminosity and measurement errors would hide the correlation. The 
strength of the correlation we can expect to measure depend on the size and quality of 
the data set. To examine the strength of the correlation we have made the following 
simple experiment. We assumed a complete positive correlation between magnification 
and residuals. The residual of a SN, the magnification of which is /i, would then be 
given by A = A M . To this residual we added normally distributed noise according to 

§ The best fit straight line is found by minimizing 
2 _ \p (Aj-a- bA^f 

x V ^ + &2 <* ' 

with respect to a and b. Uncertainties in residuals and magnifications are denoted by <7Ai and ta^i, 
respectively. 
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Figure 1. Residuals vs magnifications. The solid line shows the best fit straight line 
to the data. The dashed line corresponds to the case of no correlation. 



the measurement error of the SN magnitude and an error in the cosmological model 
used to predict the average distance modulus which was subtracted from the observed 
value. The error in the cosmological model corresponds to 10% uncertainty in fl M . This 
noise hides the otherwise complete correlation. We repeated this experiment several 
times and the result is shown as the dashed curve in Figure El The most likely linear 
correlation coefficient is r = 0.18. As can be seen from the figure a large range of 
correlation coefficients can be expected for a sample of the size considered here. The 
measured correlation coefficient (r = 0.29), indicated by the vertical solid line, is in 
accordance with the expected strength of the correlation. 

We can also asses the probability of measuring the value of the correlation coefficient 
we obtain, if there were no correlation. The solid curve in Figure [2] shows the result of 
computing the correlation coefficient for a large number of simulated data sets where SN 
residuals and magnifications are uncorrelated. The probability of obtaining a correlation 
factor of r = 0.29 or higher is 9.2%. Since the probability of no correlation is 9.2%, 
the confidence level of our detection of the correlation is 90.8%. In the case where the 
magnification was computed using SIS-profiles the probability of obtaining a correlation 
factor higher or equal to the measured value (r = 0.34) is 5.8%. 

If we assume the uncertainty in our "concordance model" to be 10% in Q-^, the 
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resulting distribution of linear correlation coefficients would be a normal distribution 
with a standard deviation of 0.01 (irrespective of the halo model used). The probability 
of finding a correlation coefficient 3 standard deviations higher (r = 0.32) or lower 
(r = 0.26) than the measured value (r = 0.29) for an uncorrelated data set is 7.2% and 
12.0%, respectively. Our results are thus not very sensitive to errors in the cosmology 
used to predict the average SN brightness. 

The evidence we find for a correlation between SN residuals and magnifications are 
weak and should be regarded as tentative at the best. However, the effect should be 
detectable at higher confidence given a larger data set. To find a correlation of r = 0.18 
(the most likely value if we model galaxies as NFW-profiles) at the 99% confidence level 
requires 220 SNe with the same distribution of redshift (0.4 < z < 1.8) and similar 
uncertainties as the small sample considered here. 

Trustworthy evidence of a correlation would indicate that our models of the matter 
inhomogeneities in the Universe are reliable and would explain some of the scatter in 
the distance-redshift relation. This scatter could then be reduced, at least partially, by 
correcting for the magnification of individual SNe as suggested by [12]. 
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Figure 2. Probability distribution functions (PDF) of linear correlation coefficient for 
simulated data sets. The solid curve shows the distribution obtained for uncorrelated 
simulated data. The dashed curve shows the distribution obtained for simulated data 
with a complete positive correlation between magnifications and residuals. The value 
of the linear correlation coefficient obtained for the data (r = 0.29) is indicated by the 
vertical solid line. The probability of obtaining a correlation coefficient higher than 
this value in the case of no correlation is 9.2%. 



